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Status of Lattice Flavor Physics 

Matthew Wingate a 

a Institute for Nuclear Theory, University of Washington, Seattle, WA 98195-1550 

This talk reviews recent lattice QCD calculations relevant for quark flavor physics. Since lattice results must be 
accurate and precise to play a definitive role in phenomenology, the focus is on unquenched results of quantities 
which can be calculated most reliably. 



1. INTRODUCTION 

In his talk which preceded mine, V. Lubicz 
did an excellent job summarizing the importance 
of lattice QCD results toward the goal of un- 
derstanding quark flavor changing interactions, 
both within the CKM framework and beyond 
the Standard Model [1]. I. Shipsey gave a sim- 
ilarly inspirational talk describing the CLEO-c 
program, which will measure many quantities 
which are within reach of precise lattice predic- 
tions [2]. That comparison will generate confi- 
dence that lattice results for experimentally inac- 
cessible quantities are correct. Those two talks 
provide ample motivation for this talk, in which 
I attempt to present a snapshot of the current ef- 
forts to contribute to the study of flavor physics. 
These are still the early days of unquenched 
lattice QCD. Unquenched simulations are being 
pushed deeper into the chiral regime, and system- 
atic uncertainties are receiving increased scrutiny. 
Consequently, this talk should be read as a status 
report of ongoing research rather than a summary 
of completed work. 

Accurate, precise results are required. Before 
one can claim to see new physics, the Standard 
Model result must be solid: uncertainties need 
to be controllable. The lattice QCD formula- 
tion provides an ab initio method for comput- 
ing leptonic and semileptonic decay properties of 
hadrons. In practice, however, lattice simulations 
require at least one additional hypothesis. Be- 
cause the light quark masses cannot be set to their 
physical values, data obtained with unphysically 
heavy quarks must be extrapolated to the phys- 
ical limit. Chiral perturbation theory gives the 



functional forms for these extrapolations, but this 
expansion cannot be applied if the quark masses 
are too large. 

If one cannot firmly use chiral perturbation the- 
ory to extrapolate simulation results, then one 
is forced to extrapolate empirically, estimating 
the subsequent uncertainty by trying many fit 
functions. This introduces a hypothesis which 
is one step removed from the desired first prin- 
ciples QCD calculation. Current results based 
on simulations using Wilson-type fermions have 
used quark masses larger than m s /2, which is on 
the border of the chiral regime. (Domain wall 
fermions simulations have been pushed to m s /A.) 
It is hard to test these fitting ansatze without do- 
ing the lighter quark mass simulations. Of course 
the chiral regime will eventually be reached and, 
with it, greater theoretical control. 

Staggered fermions are computationally inex- 
pensive, so lighter quark masses can be used in 
simulations. Light sea quark masses as small as 
m s /10 have been used to date. This allows one to 
simulate in a region of parameter space where chi- 
ral perturbation theory converges quickly. How- 
ever, there is a price to be paid. An algorithmic 
trick is needed in order to have 3, not 4, active 
sea quark flavors. This "fourth-root" trick has 
not been theoretically justified, so one must hy- 
pothesize that staggered lattice QCD in the con- 
tinuum limit is QCD. This hypothesis is testable, 
and has passed a first round. Fig. 1 shows re- 
sults for several "golden," or cleanly-computable, 
quantities divided by the experimentally deter- 
mined value [3]. Calculations on the left-hand 
side were done in the quenched approximation, 
and those on the right-hand side used 2 + 1 fla- 
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Figure 1. Comparison of quenched and un- 
quenched results to experiment [3]. 



vors of dynamical staggered quarks (the MILC 
configurations, e.g. [4]). 

I believe that these two hypotheses are comple- 
mentary and both should be pursued until neither 
is necessary. 

Due to time and space constraints, some choice 
had to be made about which topics to cover and 
which to omit. I have given priority to quanti- 
ties most cleanly extracted from simulations and 
to results which include sea quark effects. I am 
unable to cover many interesting and important 
topics, most notably K — > t:tt matrix elements 
(see e.g. [5]) and theoretical methods aimed at 
reducing systematic uncertainties in future calcu- 
lations (see e.g. [6]). 

2. DECAY CONSTANTS 

The decay constant, /m, for a meson M pa- 
rameterizes the leptonic decay matrix element of 
the weak axial current: 



(0|ygM(p)) = f M P„. 



(1) 



Matrix elements of this form are straightforward 
to compute in simulations. An additional renor- 
malization calculation is needed to match the lat- 



tice current (s) to the continuum current. The re- 
sults presented below rely on perturbative calcu- 
lations of the matching coefficients, and in some 
cases the uncertainty due to truncating the per- 
turbative expansion is a leading source of uncer- 
tainty. 

The B s and D s decay constants require no 
light quark mass extrapolation, so they are more 
straightforward to compute than fg and fo- 
These quantities have been computed using the 
2 + 1 flavor MILC configurations. Using NRQCD 
for the b quark we found [7] 



f B , = 260 ± 7 ± 28 MeV . 



(2) 



In this work the statistical uncertainty is quoted 
first and the total systematic error second, un- 
less noted otherwise. The systematic error in (2) 
is mostly due to assuming the truncated 0{a 2 s ) 
terms in the perturbative matching contribute 
with coefficient 0(1). A calculation of fg a on 
the same configurations using relativistic heavy 
quarks is underway by FNAL/MILC/HPQCD. 

Several calculations of fs 3 have been done with 
rif = 2. Although the strange sea quark is 
quenched, one argues that strange quark loop ef- 
fects are small corrections. This may be true, but 
it takes us another step away from the ab initio 
calculation we desire. Results are summarized in 
Fig. 2. The only notable difference is between the 
JLQCD result [11] and (2). Before one attributes 
this to a difference between 2 and 3 flavor calcula- 
tions in general, I think further investigation into 
the dependence on light sea quark mass is neces- 
sary. The lesson of Fig. 1 is that light sea quarks 
are needed for fx, f n , and T splittings, for exam- 
ple, to agree with experiment. JLQCD find that 
if the lattice spacing is set using then m p and 
ro agree with experiment and phenomenology, re- 
spectively. However, it is also important that the 
T spectrum be reproduced. Comparing lattice 
spacings set with different quantities, CP-PACS 
[10] found a p /a-f « 1.15 on a coarse lattice and 
a p /ar ~ 1-20 on a fine lattice. The JLQCD lat- 
tices are similar in quark mass and lattice spac- 
ing, so I speculate that one would find similar 
results for a p /a-f (and thus for af K /ax) on the 
JLQCD lattices. If I include an uncertainty due 
to lattice spacing ambiguity in the same manner 



3 



+ MILC 2002 (Stag, m saa >0.7m B ) I ♦ 1 

A CP-PACS 2001 (ImpW, m sea > mJ2) I ▲ — 

A CP-PACS 2001 (ImpW, m sea > m s /2) I A- 

O JLQCD 2003 (ImpW, m soa >0.7m s ) 1 lQ.^.-.j 

□ HPQCD 2004 (ImpS, m >m s /4,n f =2+1) I B- 



50 



100 



150 200 
f Ss (MeV) 



250 



300 



350 



Figure 2. Summary of unquenched calculations 
of /b, (bold points) [8-11,7]. Solid points used 
the Fcrmilab heavy quark formulation and open 
points used NRQCD. All points above the dashed 
line have nj = 2. The sea quark action and 
masses are noted in parentheses. The dashed 
point below the JLQCD result is discussed in the 
text. 



as [10], then I obtain the dashed point in Fig. 2. 
Therefore, I conclude that difference between 2 
and 2+1 flavor results for fg e cannot be resolved 
yet. On the other hand, there is no advantage 
to averaging them together, so I quote the 2+1 
flavor result as the best current estimate. 

Lattice NRQCD breaks down as the heavy 
quark mass is lowered to near the inverse lat- 
tice spacing. Consequently extrapolations must 
be performed to the charm sector, incurring a 
large systematic uncertainty. It is preferable to 
simulate directly at the charm quark mass using 
relativistic fermions. This can be accomplished 
on very fine lattices, treating the charm quark 
just like the strange and light quarks, or on lat- 
tices with any reasonable spacing using the Fermi- 
lab formulation [12]. A preliminary calculation of 
fu s on the MILC lattices using the latter method 
was presented here [13]. Fig. 3 shows fD s yJmD B 
as a function of sea quark mass. They find 



f Ds = 263 



+5 
-9 



± 24 MeV . 



(3) 



The quoted systematic uncertainty is mostly due 
to higher order in the heavy quark expansion. 
This calculation computes part of the matching 
nonperturbatively [14]; the remaining 1-loop cor- 
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Figure 3. Unquenched calculation (n/ = 2 + 1) 
of fo a y/TiDs vs. sea quark mass [13]. (The bare 
strange quark mass is 0.04 in lattice units.) 



rection is small, so the quoted O{o 1 s ) uncertainty 
is negligible at this level. 

The D decay constant has also been com- 
puted in [13]. Computations were done at many 
points in the partially quenched parameter space, 
allowing a global fit to PQ^PT including fi- 
nite lattice spacing effects [15]. Fig. 4 shows 
fDy/mD/(fD s y , mD B ) along the fully unquenched 
line of partially quenched parameter space; other 
plots with m sea held fixed were shown too [13]. 
The extrapolated result obtained is 



= 1.20 + 0.06 + 0.06. 



(4) 



The error due to truncating the heavy quark ex- 
pansion cancels in the ratio, so the leading sys- 
tematic uncertainty is from the chiral fits. Com- 
bining (3) and (4), they find 



fn = 224 



+10 

-14 



± 21 MeV 



(5) 



is 



where the dominant systematic uncertainty 
from the heavy quark expansion. 

Work by the FNAL/MILC collaboration on the 
B decay constant using relativistic fermions is un- 
derway. Since presenting some NRQCD results 
last year [16], our efforts have been devoted to 
improving the signal-to-noise ratio for Im- 
pressive improvement was seen by using smeared 
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Figure 4. Unquenched calculation (n/ = 2 + 
1) of fDy/m5/{f D .y/m53 [13]. The global fit 
(dotted curve) to the data points, and to par- 
tially quenched data (not shown), includes taste- 
breaking 0(a 2 ) effects. The solid curve is ob- 
tained by subtracting these effects and yields the 
result in the text. 



sources and sinks for the heavy quark [17]. Fig. 5 
shows the ratio £ = /s sA /mB s /(/s^/ms) plotted 
vs. valence light quark mass. Current efforts will 
now be directed toward obtaining data at enough 
masses to permit a PQxPT fit as was presented 
for fjj. There are existing unquenched calcula- 
tions which quote results for fg based on extrap- 
olations using m sca > m s /2 (e.g. [11]). However, 
due to the importance of this quantity for CKM 
constraints and to the possibility of previously un- 
foreseen large curvature at small masses [18,19], 
1 prefer to wait until the chiral extrapolation is 
under firmer control before quoting a result and 
uncertainty. 

The simulations for heavy-light decay constants 
using the MILC lattices have been done using 1 
lattice spacing and 1 volume so far. Finite lattice 
spacing errors have been estimated by assuming 
terms neglected in the low energy Symanzik effec- 
tive theory appear with 0(1) coefficients. Simu- 
lations with finer lattices are underway to check 
this estimate. Finite volume effects were recently 
studied in heavy meson chiral perturbation the- 
ory [20] and estimated to be small compared 
to the other current uncertainties quoted above. 
However, they could contribute when chiral ex- 
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Figure 5. £ vs. valence light quark mass in 
units of physical strange quark mass [17]. The 
open symbols used point sources and sinks while 
the solid symbols utilized smearing for the heavy 
quark propagator. Different symbols correspond 
to different sea quark masses. 



trapolations of £ are further along. 

The predictions from heavy meson chiral per- 
turbation theory are crucial ingredients for ex- 
trapolating lattice data to the physical light 
quark masses. At small quark masses, chiral log- 
arithms become important. For example, in con- 
tinuum heavy meson xPT £ has a term which 
goes like (1 + 3g 2 )m 2 ln(m^/A 2 ) times a known 
constant, where g is the B — B* — ir coupling. 
Lattice data for this ratio will not be able to pre- 
cisely computed g because simulations are done at 
masses where the logarithm is small. However, a 
large uncertainty for g will eventually contribute 
a significant error to £. Therefore, it is impor- 
tant to determine g directly. First results, in the 
quenched approximation, were published in [21] 
(for the charm sector) and [22] (for the bottom 
sector). It is important to repeat these calcula- 
tions with 3 dynamical flavors, since only then 
are the low energy constants of chiral perturba- 
tion theory physical. 
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The f n and f K decay constants have been com- 
puted very precisely on the MILC configurations 
[23,24]. As pointed out in [25], their result for the 
ratio 



1.210 ±0.004 ±0.013 



(6) 



can be combined with experimental measure- 
ments of the decay rates for n and K to \£U^ to 
yield 



IK. 



0.2219 ±0.0026 



(7) 



where the error is dominated by lattice uncertain- 
ties. This result is in agreement with, and has 
the same uncertainty as, the 2004 PDG average 
|K S | = 0.2200 ± 0.0026 obtained from scmilep- 
tonic K decay [26]. 

3. NEUTRAL B MIXING 

Measurements of the mass and width differ- 
ences in neutral B and B s mixing can give us in- 
sight into flavor changing physics. Accurate com- 
putations of hadronic matrix elements of the rele- 
vant 4-quark operators are necessary ingredients 
in order to extract the high energy parameters. 
In the Standard Model, one needs to compute 

(W\(bd) V - A (bd) v _ A \B°) = |/Im|B B (8) 

to obtain \Vtd\ from Arrid- Note that the ma- 
trix element is parameterized by the factor 
in a normalization such that Bb = 1 in the vac- 
uum saturation approximation. The correspond- 
ing calculation for the strange B® , combined with 
the experimental lower bound for Am s , yields 
and even tighter constraint on |Vt<j|- 

JLQCD results [11] are still the state-of-the- 
art. They quote, in the MS scheme at /i = mi, 



B 7 



B, 



0.836(27) stat 



0.850(22) stat 



+0 
-27 
+ 18 

-0 



(56) 



EFT 



(57) 



EFT 



+5 
-0 



(9) 



(10) 



One expects the B factors to be less sensitive to 
chiral logarithms since their coefficient is much 
smaller than for the decay constant. Neverthe- 
less a calculation using light staggered sea quark 
masses with rif = 2 + 1 is underway [17]. 



The B s width difference is dominated by loop 
diagrams with a charm quark, so a measurement 
of Ar Bs would determine |V^,K a | 2 - For this cal- 
culation another matrix element is needed in ad- 
dition to Bb s , the scalar bag factor: 



(B°\(bs)s-p(bs)s-p\B») = - 



5f 2 m%B Ss 
3(m b + m s ) 2 



.(11) 



The best result to date is B Ss (m b ) = 0.86(3)(7) 
from [18] with nj = 2 improved Wilson quarks. 
This computation is being done on the MILC lat- 
tices, but the matching calculation needs to be 
completed. 

Even for physics beyond the Standard Model, 
lattice calculations of hadronic matrix elements 
are needed to connect the experimental mea- 
surements to theoretical parameters. For exam- 
ple, the complete set of five 4-quark operators 
for B° — B° mixing have been computed in the 
quenched approximation [27]. 

4. SEMILEPTONIC DECAYS 

Below we discuss how lattice calculations of the 
form factors governing semileptonic pseudoscalar 
states can be used to compute CKM matrix el- 
ements. The form factors parameterize the rele- 
vant matrix element of the weak current; e.g. for 
B -> tt£u, 



Mi 2 ) 



"<7m + 



Ptt+PB - 



(12) 



where q is the momentum transferred to the lep- 
ton pair. This matrix element can be computed 
on the lattice for a range of q 2 ; one cannot 
presently go to small q 2 without introducing lat- 
tice artifacts as \p„\ — > 1/a. 

dT 



1 



I Kb | 2 dq? 



247r 3 ^ 



(13) 



can be integrated over q 2 using lattice data for 
the right-hand side. Then the partial width, de- 
termined by experiment, can be inserted, yielding 
\V ub \. 

A preliminary result using NRQCD on MILC 
configurations (with m soa = m s /'i) was presented 
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Figure 6. B — > 7r£^ form factors on the MILC 
configurations with NRQCD [28]. 



here [28]. Fig. 6 shows /+ and fo vs. q 2 after 
chiral extrapolation of the valence quark mass. 
The data are well-fit by an ansatz which satisfies 
the kinematic constraint /+(0) = /o(0); the heavy 
quark scaling laws /+(<?^ ax ) oc yii, /o(<2ma X ) oc 

3/2. 



l/V m B anc ^ /+,o(0) oc m B ; and includes the B* 
pole and an effective second pole [29] . Integrating 
over the entire range of q 2 and taking branching 
ratios from CLEO [30] we obtain 

\V ab \ = (3.86 ± 0.32 oxpt ± 0.58i att ) x 1CT 3 . 



(14) 



CLEO has binned their data, so we can per- 
form the integration over the restricted range 
16 GeV 2 < q 2 < q 2 ^ to find 

\V ub \ = (3.52 ± 0.73 oxpt ± 0.44 latt ) x 10~ 3 . (15) 

Note that this last step decreases the lattice 
uncertainty but more than doubles the experi- 
mental error. The moving NRQCD formulation 
[31,32] promises to permit simulations in an iner- 
tial frame moving with respect to the B, so that 
one can satisfy Ip^l <C 1/a even with small q 2 . 
Another preliminary result, 

\V ub \ = (3.0 ± 0.6 expt ± 0.4 latt ) x lO- 3 , (16) 

was reported here [33] . This is a fully unquenched 
calculation on the MILC lattices with the Fcrmi- 
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Figure 7. D — > irli/ and D — > Klv form factors 
on the MILC configurations with the Fermilab 
formulation [34]. 



lab heavy quark formulation, where the integra- 
tion uses the same restricted range of q 2 . One 
finds nice agreement between this result and the 
partially quenched NRQCD results above. As the 
authors of [28] and [33] extend their results on 
the set of MILC configurations, results will be 
highly correlated and a joint analysis would be 
the best way to extract a single result. At present 
the results (15) and (16) share configurations at 
only one value of (m va i, wisca)- For the purpose of 
quoting a single number here, I average the cen- 
tral values in (15) and (16), but just carry over 
the same experimental and lattice errors 



\V ub 



(3.3 ± 0.7 oxpt ± 0.4iatt) x 10" 



(17) 



The larger experimental error does not let us off 
the hook; extending the lattice results more accu- 
rately to lower q 2 will reduce both uncertainties, 
as (14) indicates. 

Finalized results for the D — > 7r£^ and D — > 
Klv form factors were presented here [33,34]. 
These used the Fermilab heavy quark formulation 
and fully unquenched light and strange quarks. 
The form factors and their q dependence are 
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Figure 8. B — ► Z?£^ form factor at zero recoil vs. 
sea quark mass (in lattice units; am s = 0.04) [33]. 



shown in Fig. 7. Note the entire range of q 2 can 
be studied in the D rest frame without the pion 
having large spatial momentum in lattice units. 
Integrating the fits yields 

\V cd \ = 0.239 (20) oxpt (10) stat (24) sys (18) 
\V CS \ = 0.969 (24) oxpt (39) stat (94) sys . (19) 

The leading lattice systematic uncertainty is due 
to the neglect of higher order terms in the heavy 
quark expansion. 

Ref. [33] also presented a preliminary un- 
quenched result for the B — > Dlv form fac- 
tor evaluated at zero recoil, J-b^d{w = 1) = 
1.074(18)(16). Fig. 8 shows the light quark mass 
dependence of the form factor. Combining this 
result with experiment gives 

\V cb \ = 0.038 ± 0.006 cxpt ± 0.001 latt ■ (20) 

The 3 results (18), (19), and (20) satisfy CKM 
unitarity: 

(\V cd \ 2 + \V CS \ 2 + \V cb \ 2 ) 1/2 = 1.00(2)(10) (21) 

where the first error is experimental and the sec- 
ond is from lattice. This is the first time such a 
test could be done using only inputs from exper- 
iment and unquenched lattice QCD. 

Progress was made recently in calculating the 
K — > irlv form factors [35]. Utilizing and ex- 
tending the double ratio method developed for 
B — > Dlv decays [36], the form factors can be 
computed more precisely than before. A plot of 
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Figure 9. Quenched K — ► ntv form factor [35]. 



f+(q 2 ) in the quenched approximation is repro- 
duced in Fig. 9. Using the value extrapolated to 
q 2 = with experimental input yields 

| | quenched = 0.2202 ± 0.0025 . (22) 

This precision is comparable to other methods, 
in particular (7), so it would be profitable to use 
this method in an unquenched calculation. 

5. NEUTRAL K MIXING 

In order for experimental measurements of 
K° — K° mixing to constrain the CKM parame- 
ters ~p and 77, a precise determination in QCD of 
the matrix element 

mOAs^lK ) = lf 2 K m 2 K B K (23) 

is necessary. There have been many quenched 
calculations of B#; a review was given recently 
in [37]. Several new results have appeared since 
that review [38-43]. 

The present situation of the quenched results, 
is murky, but perhaps not dire. Fig. 10 shows sev- 
eral different quenched results for Bjf plotted vs. 
lattice K mass. The cross is a staggered fermion 
result in the continuum limit [44] . The rest of the 
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Figure 10. Quenched Bk as a function of m 2 K . 
Selected results, in legend order: [44,45,39,40]. 
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Figure 11. Bare Bk vs. mass of valence quark 
computed on the n./ = 2 RBC lattices with do- 
main wall fcrmions [47]. The physical result is 
obtained by interpolating to am va \ — 0.022, cor- 
responding to m s /2. The data at am va \ = 0.01 
are discarded due to poor plateaux. Curves are 
from a fit to NLO PQxPT. 



Table 1 

Selected quenched Bk in the continuum limit (an 
asterisk indicates a 2 point extrapolation). 



Rcf. 


fcrmions 


B K (MS,2 GeV) 


[44] 

[46]* 

[39] 


staggered 

DWF 
twisted m 


0.628(42) 
0.575(20) 
0.592(16) 


[40]* 


DWF 


0.570(20) 



points are at a single lattice spacing, but are in- 
structive nevertheless. One observes a downward 
curvature in as one decreases the K mass, 
as expected from chiral perturbation theory. Re- 
ports of linear-only mass dependence have been 
made based on data which do not disagree with 
this plot: either the mass range covered is too 
narrow or m 2 K is not small enough to resolve the 
curvature within statistical errors. 

There is noticeable lattice spacing dependence 
in the RBC domain wall fermion results indi- 
cating an upward trend in the continuum limit. 
Indeed before comparing results too closely, dis- 
cretization effects must be removed. Table 1 lists 



several results for Bk in the continuum limit (the 
DWF extrapolations both based on 2 lattice spac- 
ings). The disagreement between results is not 
too severe. A clearer picture may emerge as stud- 
ies with improved staggered fermions [41,42] and 
overlap fermions [43] are pursued further. 

More importantly for phenomenology, un- 
quenched calculations for Bk are underway. Fur- 
thest along are the nt = 2 domain wall fermion 
results [47]. Fig. 11 shows Bk as a function of va- 
lence quark mass; different sea quark masses are 
denoted by different symbols. At am va i = 0.022, 
where physical results are quoted, the data tend 
to decrease with the sea quark mass. They quote 
B K (MS,2 GeV) = 0.509(18) using only degen- 
erate mass mesons. If they use nondegenerate 
valence quarks, and extrapolate the light mass 
to the physical limit they find 0.496(17). A cor- 
related fit resolves the difference between these, 
hinting that SU(3) breaking effects are becom- 
ing measurable. These results are significantly 
lower than the quenched results in Table 1, and 
effect larger than anticipated from QxPT [48]. 
The authors caution that simulations at finer lat- 
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tice spacing (as well as finite volume and lighter 
sea quark mass) are necessary to estimate all the 
systematic effects. 

Unquenched data for Bk using improved Wil- 
son fermions recently appeared [49]. The au- 
thors note that the errors are too big and the sea 
quark masses are too large for reliable extrapola- 
tion to the physical masses. However, they argue 
that their data indicate Bk should decrease in 
the limit of physical sea quark masses based on 
the observation of the partially quenched points 
(shown in Fig. 12) decreasing at any given value 
of m va i as the sea quark mass is decreased. Given 
that the lightest ?7i sca ss 1.4m s , I do not think the 
trend seen in Fig. 12 necessarily continues all the 
way to and through the chiral regime. 

Bk is also being calculated on the nt = 2 + 1 
MILC configurations; preliminary bare results 
were presented here [42]. There is some concern 
that, after the renormalization constant is com- 
puted, this result will agree with the quenched 
JLQCD result [44] and consequently disagree 
sharply with the n f = 2 RBC result [47]. Such 
comparison, however, is unjustified until the 
matching is done for [42] and systematic uncer- 
tainties are better estimated. 

6. OUTLOOK 

Now that unquenched lattice calculations are 
able to push to lighter quark masses, well in- 
side the chiral regime in the case of staggered 
fermions, predictions for phenomenology are be- 
coming more accurate. Once the ambiguities and 
uncontrolled errors of the quenching are behind 
us, we can concentrate on more reliable estimates 
of the controllable uncertainties. 

Already unquenched lattice results for many 
important quantities are being obtained. Fig. 1 
shows the comparison for many quantities which 
are well known experimentally. Table 2 summa- 
rizes the results presented here. The precision 
of the lattice calculations for \V US \ and \V u b\ are 
comparable in precision to the Review of Particle 
Physics numbers [26], and the heavy- light decay 
constants are predictions. I believe the system- 
atic uncertainties in calculations of Jb and T$k 
will soon reach the level where they too can be 
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entered in the table. These latter 2 are especially 
crucial to constraints on the CKM parameters p 
and 77. 
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